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ABSTRACT 

We investigate the spectro-photometric variability of quasars due to lensing of small mass 
substructure (from several tens to several hundreds solar masses). The aim of this paper is 
to explore the milli/microlensing influence on the flux anomaly observed between images 
of a lensed quasar in different spectral bands and possible influence of small mass structure 
lensing of non-macrolensed quasars. We find that spectro-photometric variability may be also 
caused by lensing of small mass diffuse structure and can produce the flux anomaly which is 
sometimes seen in different images of a lensed quasar. Additionally, we found that the lensing 
by small mass diffuse structure may produce significant changes in photo-center position of 
a quasar, and sometimes can split or deviate images of one source that can be detected as 
separate from the scale from 0.1 to several milliarcseconds. This can be measured ith Gaia- 
like space astrometric missions. We point out a special case where a low redshifted deflector 
(zd ~ 0.01) is lensing a high redshifted source, for which the variability in the flux and photo- 
center (several milliarcseconds) may be detected on a relatively short time scale. 

Key words: gravitational lensing: micro, galaxies: active. 



1 INTRODUCTION 

The light from a distant quasar (QSO), along its path toward the 
Earth, can be perturbed by compact objects, as e.g. galaxies, stars in 
galaxies, stellar clusters, intermediate mass compact objects (IM- 
COs with 10 2 ~ 4 M G ) and cold dark matter (CDM) structure. This 
can cause magnification in the luminosity of a QSO and in its pho- 
tometric position, the so called lensing effects. The gravitational 
lensing can be strong or weak, depending on the position of a 
quasar and object (in the line of sight of an observer) which acts 
as a lens, if the line of sight from the observer to the source lie 
very close to the center of a massive object (galaxy cluster, galaxy, 
CDM, stellar cusp, etc.) then we expect a strong lensing effect. 
The strong lensing of QSOs can cause a significant wavelength 
dependent am plification and affect the photocentric position spec- 
trum (see e.g. IOgurill2005t ICongdon & Keetonll2005l : lTreull2010l : 
lErickcek & Law 201 lb . Depending on the image angular separa- 
tion, the strong lensing can be divided into several categories from 
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which three are the most used (see Zackrisson & Rieh ml20 1 Ot iTrevj 
2010): macrolensing (> 0.1 arcsec), millilensing (10 3 arcsec) and 
microlensing (10~ 6 arcsec). Usually, millilensing and microlens- 
ing are considered to be present in the images of a macrolensed 
QSO, however the effect of micro/millilensing may be present 
in QSOs even if they are not macrolensed, i.e. the line of sight 
from the observer to source does not lie very close to the center 
of a massive galaxy, but still a compact massive obj ect from the 
galaxy and stars can affect the Q SO light (see e.g. Izaharov et all 
|2004| ; IZackrisson & Riehml 120071) . There is a n umber of QSOs 
obser ved through a galaxy stellar disc (see e.g. iMeusinger et al.l 
2010). The possibility that QSOs are micro/millilensed se ems to 
be higher than macrolensed. As e.g. Zah arov et a 3 d2004 found 
that the optical depth for gravitational micro/millilensing caused 
by cosmologically distributed deflectors could be significant and 
could reach from 0.01 to 0.1 for QSOs with z > 2. Therefore, mi- 
cro/millilensing effect may be often present not only in the case 
of lensed QSOs, but in a number of non-lensed objects. It is well 
known that micro/millilensing can be used for investig ation of the 
structure of QSOs as well as lensing objects (see e.g. lTreull2O10l : 
Zack risson & Riehml I20T0I : ISchmidt & WambsganssfeOld for re- 
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view), and especially the y present a tool for detection of dark 
matter subhalos (see, e.g.llnoue & Chiball2005l: I chen et al] |2007l ; 
IZackrisson & Riehrrjl2007l . l20ldlErickcek & Lawl201 ll etc.). 

Mainly two effects of micro/millilensing are used to inves- 
tigate structure of quasars and lensed objects: (i) a magnifica- 
tion that may be wav elength dependent due to complex struc- 
ture of QSOs (see e.g. IPopovic & Chartas 20051 ; Ijovanovic et al.l 
120081: ISluse et al J 120121 : IStalevski et alj|2012h . and in a such way 
producing a continuum (lin e) flux anomaly in the images of a 
macrolensed QSO (see e.g. iBlackburne et all 120061 : IPoolev et al.l 
l2006l : lKratzer et alj|201 ll : |Poolev et alj|2012l etc.): and (ii) a slight 
wobbling of the observed source position caused by microlens- 
ing effect, called astrometric microlensing (photocenter variations), 
that presents a photo centroid displacement of the sou rce image 
during the microlens i ng event in the lens plane ( see e.g . Hog et al. 
1995 1 : IWalkerl [l995l; jMivamoto & Yoshiil fl995l; iDominik & Sahvj 
2000; Lee et alfeoiol etc.). In the above papers it has been shown 



that for a single lens and a single source configuration, centroid 
motion during the microlensing event follows, more or less, an el- 
liptical dependence on the impact parameter of the source and lens. 

Astrometric microlensing has been explored by a number 
of authors covering different regimes. Most of them discussed 
cases where lenses are stars in the Milky Way or dark matter 
objects in Galactic halo, acting on background stars either in 
the Magellanic Clouds or in nearby galaxies , as e.g. M31 (see 
e.g. lMiralda-Escudelll996l: iMao & Wittlll998l ; |Boden et all 1 19981 ; 
iaUl" 



Goldberg & Wozniakl 119981; IPaczvnskil Il998l: [Han & Kid 119991: 
Safizadeh et all Il999l; IDominik & Sahul |200(J; iDelplancke etal] 
200 l l ; iBelokuroy & Evanl |2002|; iDalal & Land [20031 : iLee et al.l 
20ld;|Proft et al.ll201ll : lYanoll20ll etc). The effects of Milky Way 



star s on background quasars were studied in several papers (see 
e.g. Hosokawa et al. ll997l ; ISazhin et alJl998l : lHonma & Kuravamal 
2002; Sazhin et al. 201 1) an d as well as the ef fect of local dark mat- 
ter subhalos (see e.g. lErickcek & Lawll201 ll) . There are also stud- 
ies that investigate astronomical lensing in intervening galaxy (see 
IWilliams & Sahalll995l; fchen et alj 120071) . However, changing in 
magnification at different wavelengths (flux anomaly due to mi- 
cro/millilensing) in combination with astrometric microlensing is 
rarely investigated. 

A new class of satellites will be able to conduct precise mea- 
surements of an object photo-center in different spectral ranges. 
For instance, the satellite Gaia is equipped with most sophisticated 
instruments with astonishing capabilities, producing the accuracy 
better then 20/uas in wide energy range 330 nm - 1050 nm. The 
primary objective of this mission is to create a three-dimensional 
model of our Galaxy, but also to be used in other branches of as- 
tronomy, such as detection of extra-solar planets, brown dwarfs, 
asteroids, exploding stars, testing the Einstein theory, etc. One of 
the main tasks of this mission is the definition of fixed celestial 
referential frame, the Gaia r ealization of the Inte rnational Celes- 
tial Reference System (ICRS; Mi gnard et alj d2002h ). To reach this 
purpose, a large number of QSOs (5000 - 10000) uniformly dis- 
tributed over the sky is required. QSOs are Active Galactic Nu- 
clei (AGNs) located at cosmological distances. This makes them 
to a perfect candidate for this purpose as they are objects with- 
out or with negligible proper motion. The QSO sample dedi- 
cated to this task should exhibit absolutely no proper motion and 
be completely free from any contaminant (non QSO objects), 
while tiny appa rent astrometric m otions due to inner structure 
flux variations (Popovic et al. 2012) or (micro)lensing are always 
possible dTrever & Wambsganssl 2004). As we noted above, mi- 
cro/millilensing may be present also in not macrolensed QSOs. In 



the paper IPopovic et al.1 | |2012|) we investigate photocentric vari- 
ability of quasars caused by variations in their inner structure, 
and here we give our investigation of the influence of the mi- 
cro/millilensing of diffuse small mass structure (like an open stellar 
cluster) on photocentric variation of QSOs that can also have con- 
sequence for Gaia measurements. 

The aim of this paper is to investigate spectro-astrometric vari- 
ability of quasars caused by micro/millilensing taking into account 
the complex emitting structure of a QSO. Additionally we con- 
sider the lens as a bulk of stars concentrated on a relatively small 
surface, as e.g. star clusters, simulating a massive diffuse lensing 
structure that contains between several tens to several hundreds of 
solar mass stars. We explore astrometric microlensing in combina- 
tion with spectral variation from the UV to IR spectra band. This 
could be of particular interest since both depends on specific th e 
properties of a lens and a source (see lKeeton & Mousta kas 2009). 

The paper is organized as following: In the next section we 
present the source and lens model; in §3 we give results of our 
simulations and discuss results, and finally in §4 we outline our 
conclusions. 



2 THE SOURCE AND LENS MODEL 

To investigate photocentric variation in combination with possible 
wavelength dependent flux amplification, one has to take into ac- 
count the complex structure of QSOs, but also consider the lens 
structure. Here we assumed that the radiation in the UV and optical 
band is mainly coming from an accretion disc. On the other side, 
one cannot expect that a single star in galaxies which are redshifted 
more than z>0.01 can produce significant astronomical microlens- 
ing, therefore we considered a massive diffuse structure containing 
a group of stars distributed in a small surface. 



2.1 Source model - stratified emission disc in QSOs 

An AGN has a complex inner structure, meaning that different parts 
have emission in different spectral bands, i.e. a wavelength depen- 
dent di mensions of different emission regions is present in AGNs 
(see e.g. IPopovic et alj|2012l) . It is widely accepted that the major- 
ity of radiation is coming from the accelerated material spiraling 
down toward the black hole in a form of an accretion disc. The ra- 
diation from the disc is mostly thermalized from outer regions R m „ 
to the center, with slight exception at inner radius, i.e. very close to 
the black hole, where Compton upscaterin g due to the high mass 
accretion could have a significant role (see lDone et al.l20"TlT) . 

To take into account the spectral stratification of the source we 
accepted a model of a standard relativistic, o ptically thick, geomet- 
rically thin, black body di s c model (see, e.g . Pringle & Rees|[l972l : 
IShakura & Sunvaev|[l973l ; iNovikov & Thornelll973l) using the ef- 
fective temperature as a function of the radius dKrolik|[l998l) : 



1/2U/4 



T oc R- } i\\ - (Rjn/R) ' ) 

where R in is the inner disc radius. At larger radius this equation 
could be reduced to T oc R~P, where in the standard model p = 3/4. 
Here we are going to investigate the microlensed radiation from the 
UV to IR (not in X-ray), therefore we adopted the thermal emission 
as the main mechanism. 

Taking into account the tempe rature gradient along the disc we 
adopted the formalism similar as in lPoindexter et al. (2008), where 
the luminosity of a small surface element at the arbitrary position in 
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Figure 1. Temperature as a function of the position of emitting material in 
the disc (in logarithmic scale) for different p parameters: = 0.5 - full line, 
P = 0.75 - dashed line, and /? = 1.0- dash-dotted line. 



the disc is proportional to the surface energ y density and magnitude 
of emitting surface jPoindexter et al . 2008) 



he 



(i) 



where dS is the surface element of the source. We replace T in the 
expression for the energy density with the distance R as shown in 
Fig. Hand compute the proportionality coefficient as 



(2) 



where To is the temperature at the distance Rq. To compute the 
spectral energy distribution (SED) for the disc configuration we in- 
tegrate over the whole disc area: 



j 

Js d , 



dL(A, R) 



(3) 



Using Eq.[3]we calculated the SED of the source (see Fig.[2}. 
As it can be seen in Fig. [2] (up), the SED depends on value of the 
parameter /3, and for fj = 0.5 (full line) the SED has a maximum at 
lower energy than other two SEDs (with [} = 0.75 and 1 , dashed and 
dash-dotted line, respectively). An increase of the source emission 
is noticeable at higher energies for higher values of the parameter 
P- 

The inclination of the disc with respect to the observer could 
be incorporated with additional cos(i) (/-inclination angle). In all 
cases we considered a face-on disc (i = 0). 

Using the disc model for the continuum emission of QSOs we 
were able to model emission in the spectral range that will be cov- 
ered by Gaia (see Fig. [2] down). Moreover, such model allowed us 
to explore, both flux anomaly and photocenter variation in different 
spectral bands during a micro/millilensing event. 

2.2 Lens model - diffuse massive substructure 

As we noted above, the light from distant QSOs may cross far 
enough from a galaxy that is not macrolensed, but still may be 
micro/millilensed by structures which belong to the galaxy. Sin- 
gle star, or randomly distributed stars in the galaxy will probably 
act as a weak lensing, that is out of the scope of this paper. 

However, the light may cross near the center of some massive 
structures which have significant surface density that can produce 
strong (micro/milli) lensing effect. As widely accepted, galaxies 
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Figure 2. The SED of a modeled source for three values of /J (notation is 
the same as in Fig. 1). Upper: For a wide spectral range. Bottom: For U, B, 
V, R filters (separated by dotted verticals) which are considered in the paper 
and also the spectral range which will be covered by Gaia is also shown. 



can contain substructures with few tens to few hundreds stars which 
are usually addressed as the open star clusters. They can be irregu- 
lar in shape with dimension of a few parsecs or less. Here we will 
discuss effects of group of stars usually referred as an open star 
cluster that contains from a few tens to few hundreds Sun like stars. 

Such substructures contain enough stars to be able to produce 
the gravitational lensing on milli/micro scale. Here we investigate 
such structure as a possible lens. For magnification map simula- 
tion we used the ray-shooting technique, and first we describe the 
technique in a few sentences. 

2.2.1 Generation of micro/millilensing magnification maps 

A distribution of stars in the lens plane generates microlensing 
magnification map in sour ce plane which could be computed by 
the ray-shooting technique IIKavser et a l. 1986; S chneider & Weissl 
ll986lll98llTrever & Wam bsganss 2 001 . This technique proposes 
to follow a ray of light from the observer point passing towards the 
lens and falling on the source plane. Light rays initiated with differ- 
ent angles close to the observer, will have different deflection an- 
gles at the lens plane and reach the source plane at different points. 
We place the grid in the source plane and count the number of light 
rays falling onto particular grid cells. In that way, the unique map is 
created representing the current star distribution in the lens plane. 
This map is called the magnification map and can be used to deter- 
mine the amount of magnification of source imposed by the lens. 
In Fig.[3]we present a sketch of the ray-shooting technique. 
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Figure 3. A schematic representation of the ray-shooting technique in 
which a beam of light rays is shot from the observer through the lens plane 
towards the source. In the source plane we sketch the source with concentric 
curved lines. 



2.2.2 Distribution of deflectors in the lens plane 

A microlens is often defined as a random distribution of sun-like 
stars with masses equal to M G , over the lens surface. This distri- 
bution is determined by the parameter k, which presents average 
surface density. In nature, stars are mostly distributed in clusters 
with different shapes. Therefore, considering the discussion at the 
beginning of j j2,2| here we used a random star distribution with the 
application of a scale parameter x, which could allow us to change 
the compactness of the star cluster. Also, we can simulate density 
perturbation in a form of star groups in order to produce galaxy sub- 
structures mentioned in chapter JT] Additionally, we have adopted 
a thin lens approximation, which assumes that all stars belonging 
to the lens are placed in one plane. This approximation is accept- 
able when the lens and source are at cosmological distances, so we 
neglect the dimensions of the lens itself. 

First, we randomly distributed a number of solar mass stars 
Niens m the circular area with angular diameter d t . In the source 
plane of a squared shape with dimension d s , we considered a grid of 
NxN pixels, where one pixel has dimension of dJN. The deflection 
angle has been calculated using the gravitational potential in each 
point of the lens plane as: 



(4) 



For the calculation of the gravitational potential we solved the 
Poisson equation W 2 if/ = 2k in the lens plane, and we used the lens 
equation for computing the coordinates of points (elements of the 
source plane grid) where light ray will fall. This process has been 
repeated a number of times for different incoming angles in order to 
produce the light rays distribution on the source plane, i.e. to gener- 
ate the magnification map. The magnification of a point on the map 
is higher if more light rays fall on it. Those light rays are then dis- 
tributed over all pixels in the source plane, and the number of rays 
per pixel is proportional to the magnification due to microlensing 
at that pixel in the source plane. 

Using a magnification map created in described way in the 
source plane, we considered a source of light of a finite dimension 
d src , that is also divided in a number of pixels. We accounted only 
for the light rays falling on the source itself or better to say rays 
emitted from the source surface. That radiation is additionally in- 
creased by the amount defined with distribution of caustics of mag- 
nification map in particular pixels. The emission of source pixels 
is determined by it's SED. Depending on the lens mass distribu- 
tion, in the scope of the ray shooting technique, there can be more 
then one pixel in the lens plane where light rays sent from the ob- 
server penetrate and fall on the particular pixel in the source plane. 
Finally, any pixel of source has a collection of pixels in the lens 



plane that depends on it's SED. This procedure allowed us to cre- 
ate an image of the source in the lens plane in considered spectral 
range. Luminosity of the lens plane pixel L l *"*(A) is higher if there 
are more light rays passing through it and mathematically presents 
a sum of SED's for all passed light rays: 



(5) 



We calculated the centroid shift of the image for a spectral 
filter as: 



^lZ„ plx x pix L^(X)dX 
f A X npix L<™(A)dA 



(6) 



where F denotes integration for a particular (U, B, V, R) spectral 
band and A is the whole energy range. 

The magnification for a particular source image is computed 
as the ratio of the luminosity for all pixels in a spectral range with 
the luminosity in the same spectral range without lens influence, as: 



m(F) : 



(7) 



The relevant length scales for microlensing is the dimension 
of the Einstein Ring Radius (ERR) in the lens plane, defined as: 



AGm D d D ds 



c 2 D s 

and its projection in the source plane is: 



D s 

ERR = — £ 



AGm D s D ds 



Da 



(8) 



(9) 



where G is the gravitational constant, c is the speed of light, m is 
the microlens mass. We adopted standard notation for cosmological 
distances to the lens D d , source D s and between them D ds . 

Time scales for microlensing event is described in detail in 
the book on gra v itational microlens ing (see ISchneideret "all 19921 : 
IZakharovl [l997: Petters et al 2001). Based on the sizes of the 
source (R src ) and caustic (r ca ustic) pattern we distinguish two cases, 
when R src > r caustic and R src < r cmlstic (presented in detail in 
Ijovanovic et al.l |2008)). Both cases could be expressed in a single 
form: 

R 



1>CTO 



(10) 



v ± (DJD d ) 

where R replaces the R src or r caustic . Here we used a simple ap- 
proach, since we considered that a microlensing event duration cor- 
responds to the time needed for crossing over the caustic network 
created by the lens. In this way the dimension of the caustic pat- 
terns in the magnification map determine the total time scale for 
particular event, and it can be computed by using the Eq. [10] with 
the R replaced by the dimension of map r m . We used already intro- 
duced comoving distances D s and given map dimensions in ERRs 
to calculate maps linear dimensions, and hence width of the caustic 
network. 

In all calculations we assumed a flat cosmological model, with 
Q. M = 0.27, Q A = 0.73 and H = 71 km s-'Mpc -1 . 

2.3 Parameters of source and lens 

Source is defined with the inner and outer radius. As we adopted the 
standard model for the disc we considered that the most of the radi- 
ation in the observed energy range, is coming fro m the disc part de- 
fined with the R in = 10 13 m and R ou , = 10 15 m, see Blackburne et al. 
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J20T ]]). For evaluating the proportionality coefficient a, given by 
Eq.[2] we took the temperature of To = 2 ■ 10 4 K, with the peak radi- 
ation around 150 nm, at the radius R = 3.15- 10 I3 m. The coefficient 
P usually has the value of = 3/4 and we kept that value constant 
throughout all simulations. The disc inclination can be changed, 
but in the most simulations we assumed a face-on disc orientation. 
The source plane dimension is the same in all our computation and 
equal to 40ERR. 

The lens has been assumed to have a circular shape, contain- 
ing N„ stars of Solar mass ranging from 40 to 240. We also assumed 
that lens and source are placed at the cosmological distances, with 
Zd = 0.5 and z s = 2.0 (standard lens). Those values are not con- 
stant in the case when we examine dependence on the source lens 
distance. We are confident, based on the discussion in the j j2.2l fhat 
such carefully chosen lens reflect good enough condition for the 
gravitational bound systems. Any more massive and densely popu- 
lated lens will act as one compact object with known influence on 
the distant source. 



3 RESULTS AND DISCUSSION 

As illustrations, we presented in Fig. |4] the lens with random dis- 
tribution of 40 and 240 stars (the distributions of stars located at 
Zd=0.5 are shown on upper panels, and proper magnification maps 
for a source at z, = 2.0 on bottom panels). As it can be seen from 
Fig. E] an increase of the star density brings that the lens is very 
similar to the point like one, therefore, we considered the maximal 
star number as 240 solar masses stars. The source images in U, B, 
V and R filters are shown in Fig. [5] As it can be seen the source im- 
ages are different in dimension for different filters, and there is an 
expected magnification in size, that is increasing for the lower en- 
ergies. Also, there is a noticeable change in the source structure due 
to the microlens effect that depends on the lensing system, as well 
as on the emission region of the source. In the case of N s = 240, 
which we considered as the massive lens, the images are composed 
from several bright similar structures in all filters, but one can see 
that the dimensions of bright structures are different. The images 
of almost point like source, in this case, are extensive and could be 
with dimensions of 0.1 mas (see Fig.[5]right). It means that e.g. a 
non strongly lensed quasar, due to lensing by diffuse small mass 
structure, may show a complex shape with different bright spots in 
the center and arc-like structures on the edge. 

In some cases, one can expect that an image of the source is 
splitted into several components. Especially, in the case with low 
redshifted lens and high redshifted source. As an example in Fig. [6] 
we presented the images in different filters in the case where source 
is located at z s = 2.0 and lens at Zd = 0.01. It is interesting that 
the lensed source can be seen as several images distributed across 
several milliarcseconds. This can affect a relatively strong jitter in 
the photocenter of a source, since random distribution of caustics 
contribute to the source amplification. 

In the next chapter we present the influence of different pa- 
rameters to the magnification and photo-center offset in the case 
where the lens is a stellar system with mass from 40 to 240 solar 
masses. 



3.1 Population of the stars in the lens 

It is well known that a single lens microlensing of a source can 
cause a centroi d motion, or wobbling of the photo-center of the 
source (see e.g. lHog et alJl995l ; IWalkeJl995l;lMivamoto & Yoshiil 
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Figure 4. Upper panels show a randomly star distribution in the lens plane, 
a) case of n\ ens = 40 (left panels) and b) case for n\ ens = 240 (right pan- 
els). In the bottom panels we give the appropriate magnification map in the 
source plane for specified n\ ens , respectively. Horizontal parallel solid lines 
present the paths of the source for which we calculated magnification light 
curves and variability of the photo-center positions. These plots are given 
for a standard lens system with zd = 0.5 and z, = 2.0. 



1 19951: 1 Dominik & Sahu 200f]: lLeeetaT]|20ld etc.). Here, we con- 
sidered an open cluster configuration of the lens, taking different 
number of stars, ranging from 40 to 240 Solar mass stars, on the 
same surface. We performed the simulations for a typical gravita- 
tional lens system, taking Zd = 0.5 and z s = 2.0. The results of 
our simulations are presented in Fig. UJ and in Table [TJ Fig. UJ (up- 
per panels) presents the variation of magnification for six different 
star numbers in the lens, ranging from 40 to 240. They are grouped 
and randomly distributed in a circular shape with diameter around 
0.08 mas, as it is shown in Fig. [4] As it can be seen from Fig. |7J 
(upper panels), a significant difference between magnification in 
different filters is present for massive lens (A'., > 100), and can be 
very prominent, i.e. the magnification in the U filter is significantly 
higher than in the R one, and as it can be seen in TableQ] the magni- 
fication can be around 1 .4 times higher in the U filter than in R one. 
This can also affect the flux anomaly between different images of a 
lensed quasar. On the other side, the light curves in different filters 
have similar shapes, and it is expected that with a higher lens mass 
(higher stellar density) the lens acts as a point-like lens, i.e. the light 
curves have a Gaussian-like shape (see Fig.|7j upper - panel right- 
down). In the case of low density star clusters we can distinguish 
separate peaks in the magnification curve induced by the absence 
of caustics over the source path. 

The photocenter variability for considered cases is given in 
Fig. I7J (bottom panels) and maximal photocenter offset in Table Q] 
As it is expected, the photocenter variability is higher for massive 
lens. Fig. |7j (bottom panels) also shows that in almost all cases a 
decrease of the photocenter offset with approaching the center of 
the diffuse mass lens is seen, and the drop is higher for the mas- 
sive lens. As it can be seen in Fig|7j(bottom panels) and Table[T]in 
difference with the magnification, there is no significant photocen- 
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Figure 5. Images of source in four different energy channels U[332-398nm], B[398-492nm], V[507-595nm] and R[589-727nm]. Lens for this case is presented 
in Figure[4] Left four panels present the image for the case of ni ens = 40, while four on the right side are for the «/ c „, = 240 stars, with mentioned z.d = 0.5 and 
z s = 2.0. 
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Figure 6. The same as in Fig. [5] but for zd = 0.01 and z, s = 2.0. 



ter offset difference among different filters with exception for the 
last two cases of 200 and 240 stars. We conclude that for cluster 
with higher star population this mutual band deviation will not in- 
crease significantly, regarding the single lens behavior of the lens. 
The maximal photo-center offset is small, around 0.06 mas for the 
case of a high density cluster (see TableQJ. 



3.1.1 Different lens and source redshifts 

It is interesting to explore a situation where star clusters, containing 
between several tens to 240 stars, are located at low redshift and act 
as lenses to different redshifted quasars. Therefore, we performed 
numerical simulations, assuming that the lens is a low redshifted 
star cluster (id = 0.01). The linear dimension of the lens is held 
constant, consequently when moved closer to the observer it's an- 
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Table 1. Data for the maximal centroid shift in [mas] and magnification per energy channel for different number of stars in the lens, ranging from 40 to 240. 
First column gives the number of micro lenses n\ ens , next four are centroid shift in U, B, V and R energy channels, and next four are magnification variability 
in same channels. Calculations are made for the standard lens and source distances. 



ft lens 


K 


U 


dfd 

B 


nas\ 

V 


R 


U 


mag 

B 


J.H 

V 


R 


40 


0.0940 


0.011 


0.011 


0.011 


0.011 


4.6 


4.4 


4.2 


3.9 


80 


0.1620 


0.026 


0.025 


0.024 


0.026 


8.5 


8.1 


7.8 


7.1 


120 


0.2182 


0.035 


0.035 


0.035 


0.035 


10.5 


9.8 


9.1 


8.2 


160 


0.2650 


0.042 


0.043 


0.042 


0.043 


12.2 


11.2 


10.4 


9.0 


200 


0.3070 


0.050 


0.049 


0.051 


0.052 


12.7 


11.8 


11.2 


9.7 


240 


0.3440 


0.058 


0.059 


0.058 


0.058 


18.6 


16.9 


15.4 


13.0 



a 




b 


n,„=80 


c 

A/ 


n,„„ = 120 

A_ 


d 

J 


n„,= 160 


e 


n,„ = 200 


/ 


1 



a 


n, t „ = 40 


b 


^ = 80 


c 


n„„ = 120 


/A 

d 




e 


= 200 


f 


V n„„ = 240 



Figure 7. Magnification (upper) and centroid shift (bottom) variation for 
different star numbers in the lens. Different line colors represent observed 
energy channels, f/[332 - 398ran] - violet line, B[398 - 492nm] - blue line, 
V[507 — 595nm] - green line and R[589 - 727nm] - red line. 
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Figure 8. Same as in the Figure|7] but for different values of source redshift 
z s and the lens contains 80 Solar mass stars. 



gle diameter is increased. The lens contains 80 Solar mass stars 
concentrated in a circular shape with diameter 2.7 mas. The dis- 
tance of a source is changed from z s = 0.05 to z s = 2.0. 

Results of simulations are presented in Fig. [8] and maximal 
photo-center offset and magnification in Table[2] As it can be seen 
from Fig. [8] and Table [2] the amplification is almost constant after 
Z s = 0.5 source redshift and photocenter offset also stays almost the 
same, taking values around 3-5 mas. Therefore, one can notice that 
a low redshifted star cluster acting as a lens will produce similar 
amplification and photocenter offset on all high redshifted quasars. 
In all cases a flux anomaly is presented and as it is expected, due to 
dimension of the UV emitting region, the highest amplification is 
seen in the U filter. 

On the other hand, the time scales of the lens passing the lens 
across the source are smaller, consequently a high variability in the 
flux, can be detected in a period of several years (see Fig. |8}. This 
indicates that high redshifted sources projected very close to the 
center of a low redshifted galaxy can be affected by gravitational 
microlensing by small mass clusters, and therefore, they are not 
suitable candidates for astrometry. 



3.2 Influence of impact parameter of the microlensing event 

All above simulations have been performed for a source crossing 
the center of a stellar cluster (lens), but there is possibility that 
the projected distance between the lens and source are different. 
Therefore we simulated cases of different impact parameters b, for 
a small mass (40 stars) and more massive (240 stars) lens, taking a 
standard lens system (zd = 0.5 and z s = 2.0). 

The parameter b is given in relative units in comparison with 
the source plane half height and has been changed from (when the 
lens and source are centered) to 0.5 (lens located at the half height 
of the lens plane). We simulated transition of the lens taking paths 
as presented in Fig. [4] (horizontal lines, the line crossing the center 
corresponds to b=0 and the outer one to b=0.5). 

For this simulations we selected two extreme cases with 40 
and 240 stars in the lens in order to have a clear picture of what 
happens during the variation of this parameter. We showed them 
in the form of Tableland Figs.l9landll0l For example, in Fig. [9] 
(upper panel) we see that for the case of lens with 40 stars, there 
are more or less similar diagrams with casual peaks of the magni- 
fication having similar height. Those are caused by a rare caustic 
distribution in the magnification map for a small number of stars. 
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Table 2. Same as in the Table[T] but for different values of source distance defined with z, s , ranging from 0.05 to 2.0. The results are given for the case when 
the lens contains 80 Solar mass stars (k = 0.162) located at close distance of zd = 0.01. 
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One can conclude that one peak in lightcurves presents a variation 
caused by a single caustic in the map (caused by a group of stars), 
with relative variation as much as 5 to 7 times. On the contrary, in 
the second case (Fig.llOt. with a more massive lens (240 stars) we 
have the proportionally higher relative magnification of the order 
of 20 times, caused by an action of a bulk of overlapping caustics 
which gradually decreases with an increase of parameter b. For the 
last two cases ( e) and f) in Fig. [TO) the magnification is neglected. 
Such behavior suggests that the lens acts as a compact object with 
the mass equal to the sum of masses of all stars in the cluster. In 
both cases relative deviations per observed energy channels is at 
maximum for the moment of the highest magnification with the 
maximum in the U channel and decreased for another channels. 
These observed properties are more clearly given in Table [3] col- 
umn magn. 

For the case of the centroid shift variable D cs presented in 
Figs.|9]and[T0](bottom panels) for both cases of 40 and 240 stars we 
can notice a lot of similarities with main difference in the amount 
of shift which is 2 to 3 times higher, for latter. For the both cases we 
have two peaks which are evolving with increasing of the parameter 
b. 

The most interesting is that in the case of low mass cluster, 
where one can expect almost same magnification and photocen- 
ter variation in all passbands form b=0 to b=0.5. Such behavior is 
caused since randomly distributed groups of stars can be homoge- 
nously distributed within the cluster. 



4 CONCLUSIONS 

In this paper we studied the microlensing effect caused by a small 
mass (diffuse) stellar cluster (up to 240 solar masses) defined in the 
92.2| on the compact QSO source as it is defined in 32.11 The case of 
such lensing may be present in the case of a gravitational lens sys- 
tem, where one of images can be affected by the lensing of a low 
mass cluster. Moreover, this can be present in the case of quasars 
which are not strongly lensed by the central mass of a galaxy, but 
they are projected close to the center of a low redshifted galaxy. In 
the most cases, we considered a typical lens system (zd = 0.5 and 
z s = 2.0) for a low mass cluster having 40 or 80 Solar mass stars 
in a circular shaped surface with diameter of around 0.08 mas. We 
explored situations with different distances of a source, taking that 
the lensing stellar cluster is low red shifted {zd = 0.01). We in- 
vestigated variations of magnification and photo-center in different 
(U,B,V,R) spectral bands, covering the spectral range that will be 
observed with the Gaia space mission. From our investigations, we 
are able to point out the following conclusions: 

i) A small mass stellar cluster (from 40 to 240 solar mass stars) 
can significantly contribute to the flux and photo-center variability, 
not only in the case of lensed quasars, but also in the case of quasars 
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Figure 9. Magnification (upper panel) and photocenter (bottom panel) vari- 
ability of a source during the overcrossing event. Lens consist of 40 Solar 
mass stars distributed uniformly in the circular area with the diameter of 
0.04mas. In panels from a) to f) impact parameter b changes from to- 
wards maximum value equal to the half of the halfheight of source map. 
Trajectories are presented as in the Figure Distances Dj and D s are de- 
fined with it's respectful redshifts with values for the standard lens (zd = 0.5 
and z s = 2.0). 



which are projected very close to the center of a low redshifted 
galaxy. The amplification is different in different spectral filters, 
ranging from 5 to 20 times; the UV spectrum is more amplified than 
the IR one. On the other side, the photo-center variability can reach 
several milliarcseconds, that can be detected with future missions. 

ii) The lensing of such cluster can significantly contribute to 
the flux anomaly observed in some lensed quasars. In this anomaly 
one can expect that the flux in the U filter will be always more 
amplified than in the R filter. This shape of amplification can be 
used to clarify the source of the anomaly. On the other hand, the 
photocenter offset is not significantly different in different filters. 
In the case of N L < 100 lensing, one can expect that the image of a 
lensed quasar or non-lensed quasar has a structure showing brighter 
spots around the center and arc-like structures on the edge of the 
projected Einstein ring of the cluster. Extremely, in the case of low 
redshifted lens < 0.01) and high redshifted source (zd > 0.5), 
there may be several features (or point like images) separated by 
several milliarcseconds. 

iii) In a difference with the point like lens, the diffuse stellar 
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Table 3. Same as in the Tablefl] but for different values of impact parameters of the source and lens, ranging from to 0.5 of map half height. Paths of source 
over the magnification map are presented with full lines in Figure [4] We give the results for two cases when lens contain 40 and 240 stars, (k = 0.094 and 
K = 0.344 respectively). 
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Figure 10. Same as in the Fig. [9] but for cluster of 240 Solar mass stars. 

cluster microlensing can produce several prominent peaks in light 
curve during the lensing event, especially in a low mass cluster 
(up to 100 solar mass stars distributed in the considered surface 
of 0.08 mas). It depends on surface stellar distribution across the 
projected surface of the cluster, and in the low mass stellar cluster, 
the caustics caused by grouping of stars can significantly affect the 
amplification (different in different spectral filters), as well as cause 
an offset the photocenter. 

Additionally, we should emphasize a special case with a low 
redshifted small mass cluster (za S 0.01) in combination with high 
redshifted quasars (z rf > 0.5). The effect of lensing in this case is 
very important, the light of quasars may be amplified in a relatively 
short period (several years), and photocenter variability can be de- 
tected in a relatively short time. Thus, quasars projected very close 
to the low redshifted galaxies might not be good candidates for the 
reference frame objects. In this case one can expect maximal off- 
set of the photocenter of several milliarcseconds, that is comparable 
with the offset of the p hotocenter caused by the changes in the inner 
quasar structure (see lPopovic et alj|2012l) . 
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